{
 "cells": [
  {
   "cell_type": "code",
   "execution_count": 1,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "text/plain": [
       "'TCHEM_INSTALL_PATH/example/runs/AtmosphericChemistry/CarbonBond05'"
      ]
     },
     "execution_count": 1,
     "metadata": {},
     "output_type": "execute_result"
    }
   ],
   "source": [
    "import matplotlib\n",
    "import matplotlib.pyplot as plt\n",
    "from mpl_toolkits.axes_grid1.inset_locator import zoomed_inset_axes\n",
    "from mpl_toolkits.axes_grid1.inset_locator import mark_inset\n",
    "from mpl_toolkits.axes_grid1 import host_subplot\n",
    "import mpl_toolkits.axisartist as AA\n",
    "import numpy as np\n",
    "import os\n",
    "os.getcwd() "
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 2,
   "metadata": {},
   "outputs": [],
   "source": [
    "font = {'size'   : 13}\n",
    "matplotlib.rc('font', **font)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 3,
   "metadata": {},
   "outputs": [],
   "source": [
    "data = np.genfromtxt(\"full_gas.dat\", dtype=str)\n",
    "Header = (data[0,:]).tolist()\n",
    "solTchem = (data[1:,:]).astype(float)"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": 64,
   "metadata": {},
   "outputs": [
    {
     "data": {
      "image/png": "iVBORw0KGgoAAAANSUhEUgAAAasAAAEcCAYAAAB53pugAAAAOXRFWHRTb2Z0d2FyZQBNYXRwbG90bGliIHZlcnNpb24zLjQuMiwgaHR0cHM6Ly9tYXRwbG90bGliLm9yZy8rg+JYAAAACXBIWXMAAAsTAAALEwEAmpwYAABXB0lEQVR4nO2deXyU1fX/3ycbBAj7FkIgAcIqe0Q2ZRMFRYtLrbsibkVt0driWpHW2loX/HZTWyvVuv3U2roCIiiguLCpKLLIvhMCBMienN8fdyZMNkhmJpmZ5Lx9Pa9nnucuz5lxmE/uveeeI6qKYRiGYYQzUaE2wDAMwzBOhomVYRiGEfaYWBmGYRhhj4mVYRiGEfaYWBmGYRhhj4mVYRiGEfaYWBmGYRhhj4mVn4jIpSKyRESyRKSwhp/VVUTeFJHDnuMzEYmtyWcahmGEEyZW/nMQ+CswvSYfIiJtgCXAV0AnoCVwK1BUk881DMMIJ0ys/ERV56nqy8CmispF5AYRWeMZCa0SkbP8fNQdwDZVnamqh1W1SFWXq2qx38YbhmFEGCZWNYCI3AjMAK4AWgD3Av8RkW5+dDcG2CAi/xORTBH5WkSuCKK5hmEYYY+JVc3wM2CWqn6lqsWq+h6wCLgUQERGi4ie4Njh01dr4HLgBaAt8AvgWREZWcvvyTAMI2SIBbINDBEZDSxQ1Rife8eAYkqvK8UAL6jqTz3OEc1O0G2Rqh709LUKOKaqJeIkIv8F1qvqr4L1PgzDMMKZmJNXMfxgK/CAqr5WUaGqFgAZVexrNVDR9KH9lWEYRr3BpgH9RESiRaQhEOe5bug5BHgCmCkiA8QRLyIjRaSnH496GhgqIpNFJEpExgBnAf8N1nsxDMMId8JerETkbhH50rOfab+IvC0ip1ShXV8R+VhEckRkp4j82iMkvnVGicgKEckVkU0icnM1TLsKyAHmAdGe1zlAZ1X9O/AI8BzOxX0bcD9Q7b1RqvoZbs3qD8AR4E/ANaq6rLp9GYZhRCphv2YlIvOAV4AvAQFmAcOA3qqaWUmbpsB6YLGnfg9gDjBTVR/z1EkF1gD/xO2XGuk5X6qqb9TgWzIMwzCqSdiLVVlEpAlwGJisqm9XUuenuJFIO1XN8dy7D/gp0FFVVUT+AFyoqmk+7f4B9FHVYTX9PgzDMIyqE4kOFgm46cuDJ6gzDFjiFSoP84DfACnAZk+d+WXazQOuEZFYjxNECZ69UzcCREVFDR44cGAg78EwDKPesWLFigxVbeNP20gUqydxHnInWrNpD+woc2+vT9lmz3lBBXVicHubdvsWqOozwDMADRs21OXLl/thumEYRv1FRLb62zaixEpEHsetLY1U1ZPFxis7vykV3K9KnXIUFVlYPsMwjNok7L0BvYjIE8BlwFhVrTAenw97cCMnX9p6zntPUqcQOHCizgsLCykuttB8hmEYtUVEiJWIPIlz3x6rqt9Xocky4HTPPigv44FdwBafOmeWaTceWF52vaoijhw5UgUzDMMwjGAQ9mIlIn8BpuBGVQdFpL3naOJT52ER+dCn2UtANjBHRE4RkQuBu4DH9bj741NARxGZLSK9ROR64Frg0arYdeDACQdfhmEYRhAJe7ECpuE8AD/EOT14jzt96iQCXb0XqnoYN0rqACwH/gI8BjzuU2czcA5wBs5h417gZ1XdY5WZWeEWL8MwDKMGCHsHC1WVKtS5toJ73+CE6ETtPgYG+WNXZWJVUFDAjh07yM3N9afbOk/Dhg3p2LEjsbGW6NgwjKoT9mIVrlQ2Dbhjxw4SEhJISUmhTHSneo+qcuDAAXbs2EFqamqozTEMI4KIhGnAsKSykVVubi6tWrUyoaoAEaFVq1Y26jQMo9qYWPnJidasTKgqxz4bwzD8wcTKD6KioszBwjAMoxYxsfKDmJgYc103DMOoRUys/CAmJsZGVoZhGLWIiZUfREdH1xmxWrZsGRMmTKBZs2Y0adKEwYMH869//aukfObMmZx5ZtlAH5XfNwzDqAlMrPygrkwDzp8/nzFjxjBs2DA2bdrEvn37mDFjBtOnT+eBBx4ItXmGYRgl2D4rP6grYnXLLbdw2WWXlRKmSy65hOzsbK6//nqmTJkSQusMwzCOY2LlBzExMezZs4eioiKio6NPWHfD9A0cXX20VuxqMqAJabPTTl4RWL9+PRs3buSpp54qV3b55ZczdepUPvjgg2CbaBiG4RcmVn4QExODqpKZmUmbNn4lvQw5+/fvByApKalcWVxcHK1bt2bfvn0AfPzxxzRv3rxUndzcXEaOHFnjdhqGYYCJlV/ExLiPbf/+/ScVq6qOdGobr907d+6kZ8+epcry8/PJyMigTZs27Nq1i1GjRrFgQemkyjNnzmTp0qW1Zq9hGPUbc7DwA28QVu/oJBJJS0ujS5cuvPTSS+XKXnnlFUSE8ePHh8AywzCM8tjIyg98R1aRiojw5z//mcmTJ5Oamsq0adOIj4/n3XffZfr06cyYMcOCzRqGETbYyMoP6oJYAUycOJEPP/yQxYsXk5KSQuvWrXnooYd49NFHeeihh0JtnmEYRgk2svKDuiJWACNHjmT+/PmVls+cObNa9w3DMGqCsB9ZicgZIvKWiOwUERWRa09Sf6anXkVHW0+d0ZWU9zxR3z7PoHnz5iXecoZhGEbNEgkjqybAGuB5z3EyHgXKbh56BVBVLasufQDfuElVHiq1adOmToysDMMwIoGwFytVfQ94D0BE5lSh/lGgZBeuiCQDpwNXVVB9n6pm+GOXiZVhGEbtEfbTgEFgKnAIeKOCsuUisltEPhSRMdXp1MTKMAyj9qjTYiUiUcB1wPOqmudTtBv4KXARcCGwDvhQRM44QV83ishyEVnu3QxsYmUYhlE7hP00YIBMBJKBf/jeVNV1OIHyskxEUoA7gcUVdaSqzwDPAKSnp2ubNm3IyMiguLiYqKg6rfmGYRghp67/yt4IfKqq31ah7udAlWMjtWnThqKiIg4dOuSvbYZhGEYVqbNiJSIdgHOBv1exyQDc9GCV8MbWs6lAwzCMmifspwFFpAnQzXMZBXQSkQFApqpuE5GHgSGqOq5M0+uAY8D/q6DP6cAW4FsgDrgSmIxbw6oSvmLVo0ePKr8fwzAMo/pEwsgqHVjlOeKBBz2vZ3nKE4Guvg1ERHBegC+qanYFfcbh9mN9DSwBRgLnqup/qmpUXRlZnSytvWEYRjhQrZGViGzz4xmKE4I1frRFVT8C5ATl11ZwT4FKo7Cq6iPAI/7Y46UuiNX8+fM5//zzufvuu3nxxReJj4/nnXfe4aabbmLTpk08+OCDoTbRMAwDqP40YEfcBt2q/kJH4abY4qr5nLCnbdu2AOzduzfElvhPVdLap6SkhM5AwzAMD/6sWc1S1S+qUlFEYqg4ckTE06BBA1q2bMnu3SfxyZg+HVavrg2TYMAAmD27SlWrmtb+hhtuCK6NhmEYflDdNav7ge1VrayqhZ42O6v5nIggMTGRPXv2hNoMv6hOWnvDMIxQU62RlapWO8mRP20ihfbt2598ZFXFkU5tU9W09oZhGOFAUL0BRaQomP2FO5E8srK09oZhRBLB3mdVqddeXcQ7slJVnLd85GBp7Q3DiCSCvc9KAURkgIi0EpFoETlXRMYG+TlhQfv27cnLy+Pw4cOhNsUvLK29YRiRwklHViLS9ETlqppVuro8BxQB+UB7XKSILBG5RFVvDsDWsCMxMRGAPXv20Lx589Aa4ycnS2tvGIYRDlRlGvAQbsRU0TyXAtFl7nVV1TMARORrVb3Q83phAHaGJe3btwdg9+7d5ZwUDMMwjOBxUrFS1epOFfr2+Uuf15EQ2qla+I6sDMMwjJqjuuGWYoDunsv1nn1UZblCRKJVtUhV53naxXI8ll+dwXdkZRiGYdQcVRYrERkIvOlzq1hELlTV1T73VFU3l22rqgVAnZsGbN68OQ0aNLCRlWEYRg1Tnam5/wOmqWqKqqYAtwB/OlkjEXncT9vCHhGhffv2JlaGYRg1THXEKkFV3/NeqOr7QJMqtBtQXaMiicTERJsGNAzDqGGqI1YFItLXeyEip+Bc1Os1NrIyDMOoearjYHEPsEhEvsG5rJ8CXFYjVkUQiYmJLF68ONRmGIZh1GmqPLJS1Q+AnsBjwONAb1X9sApNA4pDJCJniMhbIrJTRFRErj1J/RRPvbLHhDL1RonIChHJFZFNIuLXhuWOHTuSmZlJTk6OP80NwzDqJgUFsHIlPPUUXHcdnHJKQN1VN+p6BvBONZ9xdTXrl6UJsAZ43nNUlQnAVz7Xmd4XIpKKSyL5T1xyyJHAX0Vkv6q+UR3jOnbsCMCOHTtIS0urTlPDMIy6gSps3AhffOGOL790QpWX58pbt4ZTT4Vvv/X7EdVxXT8D+APQ1dNOcK7qLU/8HrTK+a8qaf8eTlgQkTnVaHpAVStbTLoZ2KWqt3mu14rIacCdQL0Sq2XLlvHggw+ybNkyioqK6NGjBz/72c+45pprQm2aYRjhyuHDTpQ+++z4kekZD8THw+DBcMstMGSIO1JSQMQdflKdkdXfcYkUl1NNxwp/hS5A/iMiDYENwBOq+rpP2TCgbEC8ecA1IhLr2RdWJZKTkwEnVpHG/PnzOf/887n77rt58cUXiY+P55133uGmm25i06ZNPPjgg6E20TCMUFNUBN99V1qY1q51oykR6N0bLrgATjvNCVOfPhAT7IQe1ROrLFX9f34+x2+h84OjuBHSJ0AhcD7wqohco6r/9tRpDywo024v7vNoDZTzRReRG4EbATp16lRy35tpNxLF6pZbbuGyyy7jgQceKLl3ySWXkJ2dzfXXX8+UKVNISUkJnYGGYdQ++/bB558fF6YvvoCjR11Zq1YwdChcdpk7n3oqNGtWK2ZVR6zeEJGrgFdVNb+azwlE6KqFZ13tMZ9by0WkNfAr4N++Vcs0lUrue/t9BngGID09vaROo0aNaNmyZcSJ1fr169m4cSNPPfVUubLLL7+cqVOn8sEHH3DDDTeEwDrDMGqF/Hz4+msnSsuWufOmTa4sOhr694err3bCNHQodOsW0FReIFRHrL4DXgTmeBINeqfyykZdr4hAhC4YfA5M8bnegxtd+dIWNxI7UN3OO3bsWKlYTZ8+ndWrV1e3S78YMGAAs2fPrlLd/fv3A8dHhr7ExcXRunVr9u3bR0ZGBueddx5xcXFkZ2fzu9/9zjIIG0akcvAgfPopfPKJO774AnJzXVliIgwbBjff7IRp8GBo1Ci09vpQHbF6ErgI/6byAhG6YDCA0lN7y4DJZeqMB5ZXZ73Ky4nEKlxp06YNADt37iyX3iQ/P5+MjAzatGlDixYtWLp0KdHR0WzcuJHLLrvMxMowIgFVN0r65BNYutSdv/vOlcXEwMCBTpiGD3fi1LFjyEZNVaE6YrVXVf3N0ue30IlIE6Cb5zIK6CQiA4BMVd0mIg8DQ1R1nKf+NUABsAooBs7DxTGc4dPtU8CtIjIbeBoYAVyLn5ucO3bsyPLlyyssq+pIp7ZJS0ujS5cuvPTSS4wbN65U2SuvvIKIMH78eKKjj/89ceTIEfr371/bphqGURXy82HVquOjpk8+gb17XVmzZk6ULr8cRoxwjhBhNGqqCtURq7dE5DbgVSDXe7NMpuDKCETo0oFFPtcPeo5/4QQmEedl6Mt9QGecMK4HrvNxrkBVN4vIOcATwE+BXcDPqrvHykvHjh3Zt28fed49BRGAiPDnP/+ZyZMnk5qayrRp04iPj+fdd99l+vTpzJgxg9TUVAA2b97MVVddxbp163juuedCbLlhGICb0lu2rPSUnjc4QWoqnHWWE6YRI5zHXlRkpxQU1Qr9CcpXFCn2ufRmDi41lSciRRVN7YnIPcAR/BO6sCM9PV19R1Jz5sxhypQp/PDDD+Tl5dGrV68QWlc9li5dyqxZs/jss88oKiqie/fu3HbbbVx33XXl6m7ZsoXRo0ezZcuWgJ65du3aiPqMDCMs2LsXFi92x8cfw5o1bqrPO6XnFaYRI9z6UxgiIitUNd2ftlUeWfmRMdiX33rOT+IjdEBtrVnVKL4bg71rQZHCyJEjmT+/8kFvXl4eDRo0AKBp06YkJCTUlmmGUb/ZseO4MC1eDN9/7+43buwE6ZJLYORI5z7euHFoba0Fgr1zq8LVuQCFLuzxitX27dsjTqxOxsqVK5kxYwbR0dEUFhaG7RqcYUQ0qrB5c2lx8rqQN20Kp5/u4uudcQYMGgSxsaG1NwQEVazquihVhjeKxbZt2xg0aFCIrQkuw4YNs6jyhhFsVGHdutLi5PUobtXKidJtt8GoUdCvn9vzVM+plliJyEJctuDvq1g/Chcp4iZV3eCHfRFB48aNadu2LZs3bw61KYZhhCtbtsDChfDhh+7szYPXrp0TpVGjnEjVAWeImqC6I6vRQHUWLcSPNhFJSkpKwI4HhmHUIfbsgUWLjguU94/Ztm1h7Fh3jBoFaWlhvb8pXPBnGvC/IlIdH+2quRtGOCkpKaxcuTLUZhiGESoOHYKPPnLitHDh8XQYzZrB6NEwfTqMG+dGTiZO1aa6YvUvP5+T4We7iCE1NZU333yTqm4FMAwjwsnJgSVLjk/rrVwJxcUuRcbpp8NVVzlxGjjQ1pyCQHWTL045ea36SUpKCgUFBRQV1XRQecMwQoKq29s0fz7Mm+ecIvLynGfe0KFw//1OnIYMAc92DyN4BD/pSD3FG+2hsLAwxJYYhhE09u+HDz5wAjV/Puz2hBjt3Rt++lMXJeKMM+rFPqdQY2IVJLx5n0ysDCOCyc93UcnnzXPi5F2HbtkSzjwTzj4bxo8Hz3YVo/YwsQoSnTt3BkysDCPi2LYN3n0X3nvPee8dO+ZCGA0bBr/5jROoQYNs3SnEmFgFiYYNG5KYmGhrVoYR7hQWugCw777rjjVr3P3UVLjmGje1N2aMixxhhA0mVkEkJSXFRlaGEY5kZMDcuU6c5s1zEctjYpzX3h//COeeCz17mkt5GGPbpIOI1yMwUhg9ejQiUi6cUrdu3ZgzZ07J9bJly5gwYQLNmjWjSZMmDB48mH/9y99dDIZRC6jCV1/BQw+5PE5t2zpX8oUL4Uc/gtdecwK2cCHceSf06mVCFeYENLLyJDq8DOgENCxTrKpaNs9UnSYtLY2ioiKKi4uJipBwKa1ateLOO+/k888/Ryr4xzp//nzOP/987r77bl588UXi4+N55513uOmmm9i0aRMPPvhgCKw2jAooLHQZcf/7X3ds3erup6fDr3/tRk+DB1soowjFb7ESkftxSRDXAKuByMk8WEOkpaUBLq1GfHx8iK2pGjfccAP/+te/ePnll7n88svLld9yyy1cdtllPPDAAyX3LrnkErKzs7n++uuZMmVKiSekYdQ62dnOa++//4W334bMTLfHafx4t+/p3HOhfftQW2kEgUBGVlOBJ1X19mAZUxEicgZwJzAY6ABMUdU5J6g/GrgdGAI0AzYCs1X1n2XqLCrfml5VDdJbEWlpaRw8eLCUWE2fO53Ve1b722W1GNB+ALMnzK5Wm8aNGzNr1izuueceLrroopLcVQDr169n48aNPPXUU+XaXX755UydOpUPPviAG264IVDTDaPqZGTAO+84gZo/30WSaN4cJk2CyZOd916TJiE20gg2gYyHWwFvB8uQE9AEN3r7OZBThfrDgW+Ai4FTgL8Bz4hI+WED9AESfY6AIsN7R1a5ubknqRleTJkyhYSEBJ588slS9/fv3w9AUlJSuTZxcXG0bt2affv2AbBixQpGjBjBGWecwdixY9nkzcVjGMFg7174299c8Nd27WDKFLcH6vrrXbijffvghRfgootMqOoogYysPgb6AwuDZEuFqOp7wHsAIjKnCvV/V+bW30RkDHAR8FKZsn2qGrS4hS1btiQqKoq8vOMzotUd6YSC6OhoHnnkES677DKmTp1act+bSHLnzp307NmzVJv8/HwyMjJK6nTo0IG5c+eSkJDAe++9xwMPPMALL7xQe2/CqHvs3Qv/+Y9zhvj4Yxd3r2dPuOceuOACF3PPnCLqDYGI1XTgPyJyACcmmWUrqGpxAP0Hk6bAjgruLxeRBsB3wG9VtaKpwWoRGxsbcSMrgIkTJzJkyBBmzZpVci8tLY0uXbrw0ksvMW7cuFL1X3nlFUSE8ePHA5CYmFhSFhcXR0yM7Yow/GDfPidQ/+//lRao++6DH/8Y+vQxgaqnBPKLst5zfq6Scg2w/6AgIpOAccAIn9u7gZ8CXwJxwFXAhyIyWlUrTIsrIjcCNwJ06tSp0ufFxMSUGllFEn/84x8ZOnQocXFxAIgIf/7zn5k8eTKpqalMmzaN+Ph43n33XaZPn86MGTNKYiJ6OXbsGHfddZe5thtVJzMTXn8dXn3VpdgoLoYePeDee51AnXKKCZQRkJjMIsxzVYnICNzU389U9QvvfVVdB6zzqbpMRFJwjhwVipWqPgM8A5Cenl7p+46JiSE/P5+ioiKiIyw8S//+/bn00ktL7bGaOHEiH374IbNmzeKRRx6hqKiI7t278+ijj3LdddeVap+fn8+Pf/xj7rvvPvr06VPL1hsRRW6uc5J48UW3UbegALp3d1N8l1xiAmWUQyIp/5KIHAVuPZE3oE/dkbjpyV+r6uwq1H8AuFRVe52sbnp6ui5fvrzCMu/93r1706hRo5N1VWcoKiriJz/5CRMmTOD6668/Yd21a9fSq9dJP2ajrlFc7Kb2/v1vN5LKynJu5ZdfDldcYWtQ9QARWaGq6f60Dfk0XU3gcXd/F5hZFaHyMAA3PRgQsbGxFBQUkJeXV6/E6rXXXmPu3LlkZGTw73//m759+/KnP/0p1GYZ4cBXXzmBevll2LkTEhLgwgvhyitdDL4Im4EwQkOgESwSgV8Ao4CWwAHgI+BxVd0TsHXuGU2Abp7LKKCTiAwAMlV1m4g8DAxR1XGe+qNxQvVX4EUR8e4ILFLV/Z4604EtwLe4Nasrgck4j8GAiImJoaCgICKdLALh0ksv5dJLLw21GUa4kJkJL70E//wnrFrl4vBNnAiPPQbnnQf16A85IzgEEsGiO7AEaAF8gtt82x63H+pqETldVQPat+QhndIbeB/0HP8CrsXtj/IN63Qt0Ai3/nSnz/2tQIrndRzwKJCE27v1LXCux00+IKKiooiLiyMnpypbwgyjDlFU5PY8/fOf8OabLjfUoEHwpz/BpZdC69ahttCIYAIZWf0ByAJOU9Ut3psi0hmY7ym/MCDrAFX9CKh0IltVr63g+tqK6vrUeQR4JFDbKiM+Pt7Eyqg/bN4Mc+a4Y9s2aNECbroJrrsOBgwIsXFGXSEQsRoD3OwrVACqulVEZuKm4eolDRs2JCsrC1WtMDisYUQ8BQUu3NFTT7nI5SIuD9Qf/wjnnw8Ny8a1NozACESs4oAjlZQd8ZTXSxo2bIiqkpeXR0P7R1uKSPI+NSpgxw545hn4xz9g927o3BlmzYJrr7VU70aNEohYrQZuE5H3fSNViBtKTPOU1zuio6OJjY0FICcnx8SqDAUFBRbdItIoLoYFC1xsvrfecrmiJk50ojVxonnzGbVCoJuC3wHWisirOLfv9sCPgTTg3MDNizyaN29OVlYWEHkBbWua4uJi9u7dS7NmzUJtilEVMjPhuefcVN/Gjc5B4le/ghtvdCngDaMW8VusVHWuJ5TRb4F7cU4QCqwAJqnq/OCYGFm0bt2aHTt2cODAAXJycjh06FCoTQorGjduTGvzCgtv1q2DJ5+Ef/3L5YsaORIefNBFNPdJIWMYtUlA8zGqOheYKyKNcC7sB1U1OyiWRShRUVF06tSJm2++mV27drF69epQm2QYJ0fVuZ0/8QS8954TpSuugJ//HPr1C7V1hhGcCBYegarXIlWWvn37smDBAvLz80sCwxpG2JGb6zbvzp4N33wDbdu6UdTNN7vXhhEmBJJ80TgBAwcOpKCggLVr14baFMMoz4EDTpQ6dYKpU53r+XPPuX1Sv/61CZURdlRLrESkSESGeF4Xe64rOwprxuTIYIBnM+SqVatCa4hh+LJjB9xxh3M5nzkTTjvNTf+tXu3cz21NyghTqjsNOIvjSQzDPkVIKElLS6NRo0a2ZmWEB+vWwSOPuNTvxcUu0vmMGS6ZoWFEANUSK1V90Of1zKBbU4eIjo6mX79+NrIyQsvy5fD737vsuw0auDBId97pRlaGEUH4vWYlIptEpH8lZaeIyCb/zaobDBw4kNWrV1vUBqP2WbYMJkyAU091G3rvuQe2bnVBZU2ojAgkEAeLFKCyCe6GQL3/FzFgwACysrLYsmVLqE0x6guffeZEavhwWLHCjaq2bYPf/tacJoyIJlBvwMqGDOnAoQD7jngGDhwIwMqVK0NsiVHn+fxzF/po2DAnUn/4g4uGPmMGNG0aausMI2Cq6w14u4hsE5FtOKF623vtc+wH/gLMrQmDI4l+/foRFxfH559/HmpTjLrKF1/AOefA0KHw5ZduJLV5swuL1KRJqK0zjKBRXW/ATcCHntfXAMuB/WXq5AHfAf8IzLTIp0GDBgwaNIjPPvss1KYYdY1vv4V774X//Q9atYKHH4ZbbzWBMuos1fUG/B/wP8Cbp2mWqm6uAbtKEJEzcBl/BwMdgCmqOuckbfoCfwaGAJnA08Bv1MfTQURGAY8DfYBdwCOq+lSw7R86dChPP/00BQUFJdHYDcNvtm6FBx6A55+HhAS3FvWzn7nXhlGH8XvNSlWn1LRQeWgCrAF+jktBf0JEpCnwAbAXOBX4GfBL4A6fOqnAe8CnwEDgYeBPInJRsI0fOnQoOTk5rFmzJthdG/WJ/fvh9tuhe3d45RW3sXfTJje6MqEy6gEBxwb0uK/3wHkAlkJVnw+0f1V9DycsiMicKjS5AmgEXKOqOcAaEekF3CEij3tGVzcDu1T1Nk+btSJyGm4E90agNvsydOhQAD777LMShwvDqDJHjrjgso8+CseOuSgTM2daokOj3hHIPqvmIvIJsBJ4GZjjOZ7zOULBMGCJR6i8zMNNIab41CmbwmQekC4iQZ2r69SpE+3bt7d1K6N6FBbCX/8KXbu6ab/x42HNGnj2WRMqo14SiOv674BWwBm4XFYXAGOBF3GOGEMCts4/2uOmAH3Z61N2ojoxQIXJlkTkRhFZLiLL9+8v61NSOSLC0KFDTayMqqHqUnT06we33AK9e7u9U2+8Ab16hdo6wwgZgYjV2TjB8v4K71DVj1T1amABbo0pVJTd/yUV3K9KneOVVZ9R1XRVTW/Tpk21jBk6dCjr16/nwIED1Wpn1DPWrHEbes89FwoK4L//hUWLXLBZw6jnBCJWicAmVS0CcgHfVd7/ELq09ns4PoLy4t26v/ckdQqBoCvK8OHDAVi6dGmwuzbqAvv2ufxR/fu7fVNPPOFc03/0I5e6wzCMgMRqD9Dc83orbh3IS7cA+g2UZcDpIuLr8DEe556+xafOmWXajQeWq2pBsA0aMmQIDRs2ZNGiRcHu2ohkcnNdJPS0NLcWdeutsHEjTJ8OlrDTMEoRiFgt5bhAvQA8ICJPi8hfgD/iHBYCRkSaiMgAERmAs7eT57qTp/xhEfnQp8lLuKzFczwBdS8E7gIe99ln9RTQUURmi0gvEbkeuBZ4NBg2l6VBgwaMGDHCxMpwqMKbb7r1qBkzYNQoNwX45JNug69hGOUIRKwe5HhIpT/iQiydC1wGvAXcVkm76pIOrPIc8Z7nrsLl0wI3HdnVW1lVD+NGSR1wETb+AjyG2wDsrbMZOAfnHLIauBf4maoG1W3dlzFjxvD111+TkZFRU48wIoENG1x4pAsvdNEmFiyAt96CHj1CbZlhhDV+7bMSkTjcKOQJAM/U2S88R1BR1Y847vxQUfm1Fdz7BidEJ+r3Y2BQgOZVmbFjxwLw0UcfcfHFF9fWY41wITvbhUR65BFo2NCNoqZNg5iAtzoaRr3Ar5GVqubj1nwCjdpeb0hPT6dx48Y2FVjfUHXx+3r3dqGRfvITl7X3Zz8zoTKMahCI2HwCDA2WIXWd2NhYTj/9dBYuXBhqU4za4ocfYNIkmDzZhUT6+GMX0699WUdUwzBORiBi9QtgqojcKiIdRSRaRKJ8j2AZWVcYO3Ys33//PTt27Ai1KUZNkpPjQiL16QOLF8Pjj8PKlXDGCWemDcM4AYEIyjc4x4Ynca7r+UCBz5EfsHV1jHPOOQeA9957L8SWGDXGO+84kXrwQedEsW6dC0BrEfcNIyACmTSfReWZgo0K6N27N507d+add97hxhtvDLU5RjDZvBl+/nN4+20XFmnhQhgzJtRWGUadwW+xUtWZQbSjXiAiTJo0ieeee46cnBzi4+NDbZIRKN6NvQ8/DNHR8Mc/OucJ29RrGEHF3JFqmUmTJvGXv/yFjxYtYmJSEixf7n7wkpNhxAjbFBpJvP8+3Habc6S45BJ47DHo2DHUVhlGncRvsRKRTcAFqvpVBWWnAG+papdAjKuLjB49mkYNG/LOlClM3LevdGFMDJx3HtxzD6Snh8ZA4+Rs2eLWof77X7eZ94MP4Myy0bsMwwgmgThYpAANKilrCHQOoO86S8NXXmF8bi5vHziA/vnPLqLB3r2wdKlb8/j4Yzj1VLj6ajh4MNTmGr7k5cFDD7k9U/Pnw+9/D19/bUJlGLVAoO7llTlYpAOHAuy77vH//h9cdx0X9OnD9qIiPh88GLp1g7Zt3RTgo4+6hfp77oGXX3ZeZR98EGqrDYB586BvX7jvPhcuae1aF9fP1qYMo1aolliJyO0isk1EtuGE6m3vtc+xHxePb+6Je6tnrF8PU6fCsGFMXrCAuLg4Xn311fL1mjZ1f71//jm0bOnyGz3yiIuEYNQ+27bBxRe7/w8Ac+fC669Dp06htcsw6hnVHVltAj70HIILFPthmeMN4HbghuCZGeGowvXXu7/CX32VZu3bM2HCBF577TWKi4srbjNokBOsiy92f8FfeqmLL2fUDvn5bpqvVy+Xufehh+Cbb+Dss0NtmWHUS6rlYKGq/wP+B84NG5jliWBunIhXXoElS+CZZ0q8xX7yk5/w1ltv8cknn3D66adX3K5xY9d28GC46y7Yvt3t4zGPwZrlgw+cl9+6dS5U0hNPQEpKqK0yjHqN32tWqjrFhKoKFBbC/ffDgAFw3XUlt8877zwaNmxY8VSgLyLwq1+5qaeVK93a1pYtNWpyvWXbNrjoIjjrLPf/7d13Xd4pEyrDCDkBOViIyDUiMldEvhORTWWOH4JlZETz+utuH87997tNox4SEhI477zzePXVV8nLyzt5Pxde6HIf7dsHw4bBqlU1aHQ9w+vl17On2zv129+6ZIie8FiGYYQev8VKRO4HnsMlOVwNfFzmWBwE+yKfRx91P4KTJ5crmjp1KhkZGbz11ltV62vkSOfiHhvrgqKap2DgvP8+nHKK8/KbONF5+d17r8s5ZRhG2BDIyGoq8KSq9lPVyz3TgqWOYBkZsaxaBStWwC23QFT5j/rMM88kOTmZf/zjH1Xvs3dvWLYMunRxf/k//3wQDa5HbNni/oA45xz3/2bePHjjDehs2wMNIxwJRKxaAW8Hy5CTISLTRGSziOSKyAoRqcQrAURkpohoJUdbT53RlZT3DJrRzz4LDRrAFVdUWBwdHc11113HBx98wNatW6veb1KSSz0xahRccw387nfm2l5Vjh1z6Tt69XIj04cfdht7zzor1JYZhnECAhGrj4H+wTLkRIjIT3CpSH4HDAQ+Bd4Xkco2uzwKJJY5PgY+UtUyMY7oU6behqAYnZ8PL73k1ppatKi02pQpbgD6z3/+s3r9N2vmXKqvuMJNW02b5pwCjIopLnaj0O7dXfqO88+H7793XpYNKgvEYhhG2KCqfh1AN+Br4GqgNU74Sh3+9l3Bsz4H/l7m3gbg4Sq2TwaKgMt97o3GbWxuXV17Bg8erCdl7lxVUP3f/05a9ZxzztF27dppbm7uyfstS3Gx6l13uWedd57qsWPV76Ous3ixanq6+4xOPVV16dJQW2QY9RJgufqpA4FEXV/vOT9XmQ4ShKjuIhIHDMaNlnyZDwyvYjdTceGf3qigbLmINAC+A36rqosqseNG4EaATlWJXvDGG9CkSZWml6ZPn85ZZ53FK6+8wjXXXHPyvksb5qaykpPd3qCxY+F//4N27arXT11k0ya3ofr1193+thdegMsvr3D90DCM4HD06FF27txZ7ti1a1dA/Yr6udYhIjM5SfJFVX3Qr85LP6cDsBMYpaqLfe7/GrhCVXucpH0UsAV4Q1Vv97nfAxgDfAnEAVcBNwOjfZ9TEenp6bp8+fLKKxQVQWKiE45XXjnxG8SNbvv27UtMTAyrVq3ybriuPv/7H1x2mQvT9J//wJAh/vUT6ezb5wT8r391kexnzIA774RGjUJtmWFELMXFxezdu7dCIfI9srKyyrVt1qwZSUlJfPfddytU1a+UEpGUfLGsMEoF9ypiIm4asJTLnaquA9b53FomIinAnQTqdr9yJezf79ZFqoCIMH36dG644QY++ugjxvibYfZHP3KegpMnw+mnw9/+Vmojcp3n8GG3VeCJJyAnB6ZMcetTSUmhtswwwhpVZf/+/Wzfvp3t27ezY8eOktfeY+fOnRSWWRePjo4mMTGRpKQkevXqxZlnnklSUlKpo0OHDjRp0gTA/z/EiYzkixm49ab2Ze63BfZWof2NwKeq+m0V6n4OXFo98yrAu/+pGqkjrrjiCu655x4efvhh/8UKoH9/l9DxJz9xgXOXL3dJAetyVuKcHPjzn10sv8xMlwhx1iyXa8ow6jmqSmZm5gmFaMeOHeTn55dqFxcXR8eOHenYsSMjR44kOTm55NorRG3btiXaJ9hBTRKQWInIQOB+4AygOTBEVVeKyO+AxaoacOR1Vc0XkRXAeOA1n6LxVLwG5WtfB+Bc4PoqPm4AsNsPM0szf74Lr9S2bZWbxMfH86tf/Ypf/vKXfPLJJ4wYMcL/57dq5aKD3323G2ksXuw8E/v187/PcOToUXjqKSfGe/a4yOgPPeSCABtGPSEvL49t27axdetWtmzZwtatW8sJUU5OTqk20dHRJCUlkZyczJAhQ7joootKxCg5OZnk5GTatGlDVBit7wayZjUSWICLxL4AuBVI94jVb4FTVHVyUIx0rusvANOAT3BrS1OBPqq6VUQexgnluDLt7gN+CSSqanaZsum4taxvcWtWVwJ3ARep6n9OZM8J16yOHnVrRrffDn/4Q7Xe57Fjx0hNTWXAgAHMnz+/Wm0rZe5cuPZal8jx9793CR7D6AvoFwcPwp/+BE8+6UZS48bBr3/tonoYRh3j2LFjpYSo7Ovdu0v/fR0VFUViYmKJ6JQVoeTkZNq1a1drIyJfRKT216yA3wPzgMlANE6svKzEubQHBVV9VURaAffh9kKtAc5RVe9O2kSgq28bcZOjU4EXywqVhzich2ESkIMTrXNV9b2AjP3sMygocM4V1aRx48bBG115mTDBpba4/nq44w547TW3ltW/VrbIBZctW+Avf4Gnn4YjR+C889wes9NOC7VlhuE3hw4dOqEYZWRklKofGxtLp06d6Ny5MxMnTqRz58507tyZlJQUOnfuTFJSErGxsSF6NzVHICOrbOBCVZ0rItFAAcdHVmcA81S1Ti6UnHBk9eCD7jh40G3crSbe0VWfPn1YuHBhQAuSpVB1rtt33ulGIz//uRuN+GFjraIKH38M//d/zttRxOX4uvvuyBRco95RWFjI9u3b+eGHH9i0aVPJ4b0+dOhQqfrx8fElwuMrQt4jMTExrKbnqkOoRla5QGW+wInA4QD6jlw++cSlP/dTBBo3bsyvf/1rbrvtNt566y1+9KMfBccuEbj6apg0yf3QP/44PPeci+Bw663h59a9f79bZ3v2WTcybNnSpUr56U8tS68Rdhw+fLhSMdq6dStFRUUldePi4khNTaVLly4MGzaM1NTUUqLUunXr4P2RWocIZGT1Fs6pwuu6VgAMVtVVIjIfyFDVy4NiZZhR6ciqqMiFVrriCjfV5icFBQX079+fgoICvv32W+Li4gKwthJWrXKRxt97z+0Ju+02uPHG0CZ2zMtza2xz5sA777jwUYMGuVBSl19etz0ajbBGVTlw4AAbNmxg/fr1bNiwgQ0bNpQIU2ZmZqn6rVu3pmvXrnTp0qXk8F536NAhJOtF4UAgI6tAxKo/ztlhC/A6zivwT7h4gYOBUz17meoclYrVV185L8AXXoArrwzoGXPnzmXixIk89thj3HHHHQH1dUKWLnXTlgsWuLQYV17pguMOH147jhhHjjjBfPNNdz5yxEXf8NrRt2/N22AYHrKyskqEyCtK3vPBgwdL6kVHR5OSkkLXrl1LiVLXrl1JTU2ladOmIXwX4UtIxMrz4EHAH3Gu69FAMbAEuENV62x2wErF6tlnnSPDunUuYGqAnHPOOXz66aesXbuWxMTEgPs7Id9+69aFXnjB7Vvq2BF+/GM4+2yXR6tx4+A8JzvbOaEsWgQLF8IXX7gRVJs2bjPzBRe4/Wl1cIHYCA8KCgr44YcfWLt2bTlB2rNnT6m6nTp1Ii0tje7du5ecu3fvTkpKSp10YqhpQiZWPgY0BFoChyrxvKtTVCpWt97qInsfOhSUUcmGDRvo27cvkyZN4vXXXw+4vypx5Ai8/bYLEzV3rvNsjI2F9HQ3auzb16XXSEqC9u1d/MOy8+v5+c6JY+tW58G3ebNLw7F6tRPy4mL3+aSnw5gxcO65biRXT6dGjJohJyeH9evX891337F27dqS84YNGygoKCip17Zt2xIR8hWmrl270ijc1nIjnFBNA84AOqrqbRWU/R+wXVX/6FfnYU6lYjV8uPth//jjoD3r4Ycf5p577uHNN99kcgXZhmuUY8ecw8jChW66cM0aF9LIlwYNIC7OxeATcfvMyuyEB1yg3YEDneANGeLCQdlUiREEsrKyWLt2bSlB+u6779i8ebM3wwJRUVF07dqVXr160atXL3r37k2vXr3o3r07zcLdI7YOESqx+h54TFX/XkHZdcCdqtrbr87DnArFqqgIEhKck8Ls2UF7VkFBAenp6WRkZPDtt9/SvHnzoPVdbVRhxw6XB2rPHti713ntFRS4Q9WNtJo2hebNndde587uMGEyAqSwsJD169fz9ddf880335ScfROXxsXF0aNHj3KilJaWRsOGDUNovQGhc13vROWJCjcB9Ss/+Lp1bq0nyKF+YmNjefbZZxk6dCg333wzL7/8cujcWkXcCCk5OTTPN+oFqsqePXtKCdLXX3/Nd999VxK/LiYmhp49ezJ8+HBuuukm+vTpQ69evUhNTSUmJhJCnhrVJZD/q9m46A8V0RHIC6DvyGP1anceODDoXaenpzNr1izuvfdeJkyYwLXXXhv0ZxhGKCguLmbjxo2sXLmy5Pjqq69KRW1ISkqib9++nHXWWfTt25d+/frRs2fPmtnSYYQtgYjVEuCXIvK6qpYIkyeR4S885fWHb7916zY1FOl7xowZzJ8/n1tvvZURI0aQlpZWI88xjJqisLCQdevWlYjSihUrWLVqFUePHgWgQYMG9O3bl8mTJ9OvXz/69u1L3759aRXKvX9G2BDoPqtPcSk8/o1LkJiECwjbChihql8Fyc6wosI1qwsucFOB331XY8/dvn07/fv3JyUlhaVLl5qnkhG2FBcX8/333/P555+zfPnykhGTN/p3fHw8AwYMYPDgwQwaNIhBgwbRu3dvcwev44RkzUpVvxKRMbhgsDOAKNw+q6W4yOV1Uqgq5bvvanwDa3JyMs8//zznn38+N954Iy+88IKFZTHCgr179/L555+XHF9++WVJxtiEhAQGDRrEzTffXCJMPXr0qLdRHAz/CGglUlW/AM4QkXigBXBQVXNO0qzukZcHGze6hIc1zKRJk/jNb37Dfffdx8CBA/nFL35R4880DF9yc3NZuXJlKXHasmUL4CI79OvXj8svv5yhQ4dy2mmn0b1794gNvGqED0Fxm/EIVP0TKS/r17uNrr1rx1P/nnvuYfXq1fzqV7+id+/eTJw4sVaea9RPDhw4wNKlS0uOFStWlGyqTU5OZujQodx6662cdtppDBo0yKanjRrBfDyDgXedqpbESkR47rnn+OGHH7j44otZtGgRQ4YMqZVnG3UbVWXLli0lwrRkyRLWrl0LuD1Mp556KnfccUfJqKnGw4AZhodA09pfA1yG23NVdsedqmrX8q3qIGvXuj1IQYgHWFWaNGnCe++9x/Dhwzn33HP55JNP6F6LzzfqBqrKt99+y8cff8ySJUtYunQpO3fuBKBZs2aMGDGCq666ipEjR3LqqafaxlojZPgtViJyP/AgLmvvamp4X5WITMOToh6X1Xe6qlboHi8iKcDmCoomqupcn3qjgMeBPsAu4BFVfaraxm3c6DbK1vI/5Pbt2zNv3jxGjBjB2WefzUcffUTnzvVrL7ZRPVSVDRs2sHDhQhYtWsSiRYvYv38/4PYznX766YwcOZLTTz+dPn36mBOEETYEMrKaCjypqrcHy5jKEJGfAE8C03DehtOA90Wkt6puO0HTCYCvV2JJ0hkRSQXeA/6Jc7cfCfxVRPar6hvVMvCHH6BraAaRaWlpvP/++4wbN47Ro0ebYBnl2LJlC4sWLWLhwoUsXLiQXbt2AU6cJkyYwJgxYxg9ejQpKSnmXWqELYGIVSvg7WAZchLuAOb4xCG8TUQmAD8F7j5BuwOquqeSspuBXT6BeNeKyGnAnUD1xaq2g8z6MHjwYBYsWMD48eMZPXo0ixYtIiUlJWT2GKFl//79fPDBB3z44YcsWrSIzZvdJEObNm0YM2YMY8eOZezYsXTr1s3EyYgYAhGrj3GJFhcGyZYKEZE4XDLHR8sUzQeGn6T5fzzpSzYAT6iqb56NYZ4+fJkHXCMisapa4FsgIjcCN4LLcVNCVpYL5hqikZWX9PT0EsEaNWoUc+fOpVevXiG1yagdCgsL+eKLL5g7dy5z585l+fLlqCrNmzdn9OjR3H777YwZM4Y+ffqYOBkRSyBiNR0nBgdw02mZZSuoanEA/XtpjUvsuLfM/b3AmZW0OYobIX0CFALnA6+KyDWq+m9PnfbAggr6jPE8c7dvgao+AzwDLoJFScEPP7hziMUK3Ajrww8/ZMKECYwYMYK3336bESNGhNosowbYvXt3iTh98MEHHDx4kKioKIYOHcqDDz7IhAkTGDRokK05GXWGQMRqvef8XCXlGmD/FfXni1Rwz1VUzQAe87m1XERaA7/ChYY6UZ8V3a8cr1h161blJjXJwIEDWbZsGRMmTGDcuHG89NJLXHjhhaE2ywiQgoICPv30U95//33mzp3LV1+5pdjExEQmT57MhAkTOPPMM2nZsmWILTWMmiEQMZlFdX7U/ScDKMKNhHxpS/nR1on4HJjic72nkj4LgQNV7jWMRlZeunTpwqeffsqkSZO46KKLmDlzJvfff79FEYgwtm3bVjJ6WrBgAUeOHCEmJoYRI0bw+9//ngkTJtCvXz+b2jPqBYHEBpwZRDtO9Jx8EVkBjAde8ykaT/UcIQZQempvGTC5TJ3xwPKy61Un5IcfoHVrl3gxjGjdujWLFi3i5ptvZubMmaxYsYIXXnjBsqKGMbm5uSxZsqREoL7zbDZPTk7msssuY+LEiYwdO5amlsjSqIcEuik4EZcOZBTQEjci+Qh4/AReeP7wOPCCiHyBW4e6GegAPOWx42FgiKqO81xfAxQAq3DBdc8DbsEF3PXyFHCriMwGngZGANfiNjlXnW3bIEw97+Lj45kzZw7p6enccccdpKen8+KLL1q0izDihx9+KJnaW7RoEdnZ2cTFxTFq1CimTp3KxIkT6dmzp42eDENV/TqA7sA+IB9YBLzsOefhptjS/O27kudNA7Z4+l8BnOFTNgfY4nN9DfAdcAzIApYDV1bQ5yhgpafPzcDNVbFl8ODBWkLv3qoXXKDhzuLFizU5OVmjo6P1N7/5jRYUFITapHrJsWPH9N1339Vbb71Vu3XrpripdO3ataveeuut+s477+jRo0dDbaZh1Ai4mSu/NCCQfFZvAqcA41V1i8/9zjiX8G9VtU6u7Jfks1KFpk1h6lSYPTvUZp2UQ4cOMW3aNF5++WVGjBjBs88+S48aShZpOFSV77//nrlz5/L++++zePFi8vLyiI+PZ8yYMUycOJEJEybQLUwcdAyjJglJPitgDG4kssX3pqpuFZGZwF8D6DsyOHwYjh51oZYigObNm/PSSy8xadIkbrnlFvr168ddd93F3XffbTHfgkhWVhYLFy4sWXvaunUrAL169WLatGlMnDiR008/3T5zw6gGgYhVHHCkkrIjnvK6zTZPpKcIESsvl19+OePGjeOOO+5g1qxZvPzyyzz55JNMmDDB1kb8oLi4mFWrVjFv3jzmzZvHp59+SmFhIU2aNOHMM8/k7rvvZsKECRYGyzACIBCxWo0Le/S++mz+FfdrN81TXrfZvt2dfSNaRAjt2rXjxRdfZMqUKUybNo1zzjmHMWPG8Ic//IFTTz011OaFPXv37mX+/Pklm3K9wWAHDhzInXfeydlnn83w4cOJi6v7f7MZRm0Q6D6rd3Ax9V7FuYW3B34MpAHnBm5emBOhIytfzjzzTNasWcPTTz/NrFmzGDJkCBdffDH33HMPAwcODLV5YUNubi7Lli0rGT2tXr0acPH2zjrrLCZMmMD48eNp165daA01jDqK3w4WAJ5gsr8FBnI8osQK4H5VnRcUC8OQEgeLe+6BP/4RcnOhDoS1ycrK4tFHH2X27NkcOXKEs846i7vuuovRo0fXu+nB/Px8vvzyy5I0Gp9++im5ubklm3LPPvtszj77bAYMGGCbrQ2jigTiYBGQWPkY0AhoARxU1eyAOwxzSsTqyith6VLYsiXUJgWVQ4cO8dRTTzF79mz27t1Lnz59uPHGG7nqqqto0aJFqM2rEQoKCli1alWJOC1ZsoTs7GxEhP79+zNmzJiSVBoJYbYB3DAihVoTKxGJwk3vbVbVNZXU6QukqGptpQ+pdUrEaswYKChwglUHyc3N5cUXX+Tpp5/myy+/pGHDhkyePJkf//jHTJw4kfj4+FCb6DcZGRksW7aMTz/9lGXLlvHFF1+Qk5MDQJ8+fUrEadSoUbRq1SrE1hpG3aA2xepqnEt6X1WtKBOvN0vvGuAGVX3ZH6PCnRKx6tULTjkFXnvt5I0inFWrVvH3v/+d1157jYyMDBo3bsykSZM455xzGDduHElJSaE2sVIOHz7M6tWrWbVqFStXruSzzz5jw4YNAMTExDBw4ECGDx/O8OHDGTVqlK07GUYAqCpZeVnsOrKLnUd2sjNrJzuP7GTXkV385dy/1JpYzQfW6fGEhZXVexLooaoT/DEq3CkRq+bN4eqr4f/+L9Qm1RqFhYV8/PHHvPbaa7z55pvs27cPcHuIxo4dy5AhQ0hPT6dHjx61np4iOzubDRs2sG7dOtatW8fq1atZvXo1mzZtKqnTvn17TjvtNIYNG8bw4cNJT0+P6BGiYdQmBUUF7Dm6p5wI+V7vzNrJsYJj5dq2aNiCg3cdrDWxygCmnGyKT0TOA55T1db+GBXupKen6/IlS6BRI3joIedoUQ8pLi7mm2++YcGCBSxYsIAlS5Zw7Jj7kjZu3Jh+/fqRlpZWciQnJ9OuXTvat29P48aNq/Ws/Px8Dh06REZGBjt37mTHjh0lx+bNm1m3bh3bvN6ZHrp168bAgQMZMGAAAwcOZODAgbRvXzbQvmEYqsqh3EPlRcgrQJ77+47tQ8sk24iNiqVDQgeSmiaRlOA5mia5ez6vG8U2qtUIFgnAwSrUO+ipW3fZ68lOUo9//KKioujfvz/9+/fnF7/4BUVFRaxbt47ly5fz5Zdf8u233/Lhhx/y/PPPl2vbuHFjmjVrRqNGjYiPj6dRo0ZER0dTWFhIYWEhRUVF5Ofnc/jwYQ4dOkR2dsV+O23btqVz586cccYZdO/enR49etCjRw+6detWbUE0jLpIUXER+47tY0fWDnZk7WDnkZ2lz1nunFOYU65tq/hWJSI0qP2g0qLkObdq1IooqXmP2OqKVQbQGTiZR0EnT926y25PtpF6LFZliY6Opnfv3vTu3Zurr7665P6xY8fYuHEju3btYu/evSVHVlYWOTk5ZGdnk52dTXFxMTExMURHRxMTE0NsbCzNmjWjefPmtGjRgubNm9OqVSuSkpLo2LEjHTp0oEGDBiF8x4YRWvIK89h1ZFd58Tmyo+T1riO7KNKiUu28o6GOTTsyKHEQ53U/r5wIJSYk0jAmfEKCVVesluIimr94knrXcnJBi2z2eDKgJCaG1o4IoHHjxiUjMMMwqsaRvCPlR0MeIfKOhvZn7y/XrnFsY5KbJZOUkMTY1LEkJSTRsWlHkpp6zglJtGncplZGQ8GkumI1G1gqIk8AM1Q137dQRGKBR4GxwMigWBiueMXKRlaGYVQDVSUjO6OUCFU0LXckv3zo1daNWpcIzpCkIRUKUdMGTevkJv5qiZWqLhORXwCPAVd4vAO3eoo74zLttgJ+oaqfBdXScGPPHhCBNm1CbYlhGGFCsRaz79g+th/ezvas7Ww/vN2Jkc+03M4jO8kvKvV3PtESTWJCIkkJSfRu05uzup5VTog6JHQIq2m52qbasQFVdbaIrATuAi4AvH6/Obgswb9X1SVBszBc2b0b2raFmICSLRuGESGoKpk5mSUiVOrseV2REDWMaVgy6hmePLzkta8QtWvcjuioyA/ZVpP49UurqouBxZ6IFl739AOqZVbxgoiITAN+CSQC3wLTKxNFERkN3A4MAZoBG4HZqvrPMnUWVdC8l6p+f1KD9uyxKUDDqENk5WVVKkLec1mPudioWJKaJpHcNJlhycNIbprsjmbHz63iW9XJabnaJqBhgSc1yL4g2VIpIvIT4Elc6pGlnvP7ItJbVbdV0GQ48A3wCC4a/NnAMyKSq6ovlanbB8j0uS6/YlkRJlaGETHkFOSccES0PWs7WXlZpdoIQmJCIp2adaJ/u/5MSptUSoSSmybTrkm7iHNUiFQiZQ7rDmCOqv7dc32bJ+L7T4G7y1ZW1d+VufU3ERkDXASUFat9qlp9N/t9+6Bnz2o3MwwjuHg3tG49vJWth7aWPnteV+Q116ZRG5KbJdOtZTfGpIwpJ0QdEjoQGx0bgndkVETYi5WIxAGDcV6GvszHjaCqSlNgRwX3l4tIA+A74LeqWtHUYHkOHAALcGoYNU6xFrP36N4TilFZz7mGMQ3p1KwTnZt1ZkCPAXRu3plOzTqViFHHph3rtbNCJBL2YoVbE4sG9pa5vxc4syodiMgkYBwwwuf2btzI7EsgDrgK+FBERnvW5Mr2cSNwI0Dn5GQ4etTEyjCCQGFxIdsPb69UjLYf3k5eUV6pNs0bNqdzs850adGF0Z1H07l5Zzo361xybtu4ra0T1TEiQay8lA1iKBXcK4eIjMBN/f1MVb8o6Ux1HbDOp+oyT8T4O4FyYqWqzwDPAKT3769s3w6t62ToQ8MIKqrK/uz9bDq4ic0HN7vzoc1sPuRebz+8vVyEhfZN2tO5WWcGJQ7igp4XlBKizs0707RB0xC9GyNURIJYZQBFQFlvhraUH22VQkRGAu8Bv1bVv1XhWZ8Dl560VmGhO9vIyjAAOJp/lM0HnQD5CpL3nF1QOrZju8btSG2RyvDk4aT2TSW1eWqJGCU3S7YpOqMcYS9WqpovIitwG459E0eNB96orJ2InAG8C8xU1dlVfNwA3PTgiTGxMuoZxVrMriO7WH9g/fER0qHjI6WyDgyNYxvTpUUXurbsyvgu40ltkUqXFl1IbZ5KSvMUGsdZkGGjeoS9WHl4HHhBRL4APgFuBjoATwGIyMPAEFUd57kejROqvwIvioh3VFakqvs9daYDW3B7tuKAK4HJOI/BE2NiZdRRMnMyWX9gfbljQ+aGUqOjaImmU7NOdGnRhR/1+JETIh9Bat2ota0ZGUElIsRKVV8VkVbAfbhNwWuAc1TVG+opEejq0+RaoBFu/elOn/tbgRTP6zich2ESLvrGt8C5qvreSQ3yipWtWRkRSHZBNhszN1YoSgdyDpTUi5ZoUluk0r1Vd8akjKF7q+6ktUqjW8tudGzakZioiPj5MOoI1Uq+aDjSk5J0+a5dkJMDDW1u3Qg/CosL2XJoS4WCtD1re6m6HRI60L1Vd7q37O7OniO1RSpx0XEhegdGXaQ2ky8aAEVFLkuwCZURQlSVPUf3lBajzPWsy1jHpoObKCguKKnbrEEzerTuwaiUUaVEqVvLbiQ0qNt5Uo26gYmVPxQW2nqVUWscyj3EhgMbSgmS9/XR/KMl9RpENyCtVRp92vbhgp4XlBol2RqSEemYWPlDYaGtVxlBJbcwlx8yfyg3Slp/YD37jh0PvxklUaQ0T6F7q+6MTB5ZSpA6Nu1okbuNOouJlT/YyMrwg6LiIrYd3lahIG09tBX12ePevkl7urfqzvndzy8lSF1adKFBTIMQvgujptBiRQsVLSp9pgh3XUFZydm3XkVlxbjrotKvy15T5LHD+7pI0eJKXle3r+LA/CNMrPzBxMqoBG+0Bl9BWndgHesPrGdj5sZSuY4S4hLo3qo7wzoO45r+15QSJYvQUDHFhcVonlKcV0xxXjGarxQXFKMFWnLUynX+cSEIlsCcPB5PGBANEiVItLjX0YJE+byOFog6/rrsdSCYWPmDiVW950jeETZkbqjQ2+5w3uGSerFRsXRr2Y3urbpzbtq5pQSpXeN2EbOOpOp+qIuyiyjOKS45inKKKM72eZ1T7K5zi0sEpTivtMCc7FrzKy+juIbfqIDEChIrRMVGlbyu9DpGkBghqkGU+8GO8fxI+559XpfUqaieT/1q9VVBWam+fMXCH4HxClRUEL6rAXRhYuUPRUXQokWorTBqmLzCPDYd3FShKO0+ejzQiSB0ataJ7q26c2W/K0sJUqdmnWplP5Kq+0EvOlpU6ig+Vv5e0dEiio4dP5cIj48Q+QqP93UgQiGx7gddGriz95AGQlSc57pRFDEtYkqXla1bQdtqicvJrgP869+oOUys/KWpTdPUBfKL8tl8cDMbMzeyIXMDGw5scOfMDWw7vI1iPf4L3aZRG7q36s6EbhNKCVLXFl2Jj4336/lFuUUUHiqk6HARhYcLS45y11mViI3PUR0xiWoYRXSTaKIaRRHdKJqo+Ch3NIoitnVsyXV0vKtT8tpTp+S1t55PHyVlvoISjL/KjXqNiZW/NGsWaguMKuLdIFsiRD6CtPXQ1lIRv5s1aEZaqzSGdRzG1f2uJq1VGmkt0+jeqjst4o+PplWV4txiJySbisg6nFW50JzgWvNPvlARnRBNdNNoYhJinMA0jiIuMY7oJtGlj8bRVb5nIwgj0jCx8hcbWYUNqsqBnAMlAVZ9o31vPriZrYe3UlhcWFK/SVwT0pqnMajFIC5OvJjU6FRSNIVOBZ1ofrQ5RVlFFO4opPDQcXHZengrPxz+oZTYaEHVhCamWQzRzdw5rm0cMWkxxDQ/fs97VHidEGPCYhiYWPmPjaxqDVXlcN5hdmTtYHPmZjbu2cim/R5hOrKZrdlbOVZ8rFSbFsUt6JjfkbScNM7IOoOOmR3psLcDHXZ1oOnuplBQ/jn7PP8BIMeFxiscce3jaNSjUeXCYkJjGDWGiZW/2MjqpKiqW6A/Vslai+f1kSNH2Jm9kx05O9iZv5NdRbvYrbvZHbWbPbF72NtgL9mxpfMhNcxvSPtD7elwsAM9D/Yk8VAi7Q+66/ZZ7Wka3/S4aDT1CEjXaGIGnkBcfK8Tom2dxTDCCBMrfwnzkZWq28dRnOfciDVPj7sTV3SugSO7OJvMhEwym2RyoMkBd05w58wmmWQkZLC/2X6OxB9xRsd7DqBlbkva57UnNTeV4bnD6RDdgaS4JJLjk0lNSKVdUrvS4tL0uNhENYqKGJdwwzCqhomVv9TAyEqLlYL9BeTvz6dgX4F7vS/feYt5RyNHyniBZRdVKkbB2JMiDYSohlFENYxC45VjCcc43OwwWU2yONz6MIcbHeZw/GEONDzAgbgDZMRmkBGTwf6o/RyTY+X6i5EY2sa1pV18O3o07sGZzc6kU8tOdGrViU7NO5HcLJmkhCSL0mAYRilMrPzgUEMCGlkVFxaT/V02R1Yc4ehXR8nZmEPuD7nkbM5B8ypetPe6Gpc6EqKJbRNbIiZRDUqfS4SmQfky4uBYzDEORR8iUzI5qAc5UHyAzKJMDhQeIDM/kwN5B9ifvZ+M7AwysjPIzMksFRLIl6YNmtK+SXsSmyTStUlXEpskuuuExJL77Zu0p1WjVkRJlN+fnWEY9RPLZ+UH8W1Fc3YXQnTVg4bmbs8l478ZZPw3g6xPsyjOdcOeqEZRxHeLJ76rOxp0bkBc2zhi28a6c5tYYprHEBV74h/4vMK8ElHxHr5CU9G1b+gfX2KjYmndqHXJ0aZxG1rHty5/z/O6VXwrv/cZGYZRfwgkn1XEiJWITAN+icsK/C0wXVWXnKB+X+DPwBAgE3ga+I36vGERGQU8DvQBdgGPqOpTJ7Wlg+jeDXtp27htpXVUley12WS8mcH+N/dzdIVL5dCoVyNaTmhJQnoCCYMTiE+LL7eQX6zFHMw5WF5sjnmuc8pcZ2dwJP9Ipba0aNiilLi0adSmRGTaNG5Tcu09mjZoams+hmEEnTqffFFEfgI8CUwDlnrO74tIb1XdVkH9psAHwGLgVKAHMAc4BjzmqZMKvAf8E7gSGAn8VUT2q+obJ7Pp2ZXPcvfpd5e6p8VK1hdZZLzpRlA563NQlNgRscT+LpbiUcUcaHOAddnrnNjs2E/G+vIjn8yczFKRE3xpFNuolOh4cxWVFRyvOLWMb2npxw3DiHgiYmQlIp8DX6vqDT73NgCvq+rdFdT/KfAHoJ2q5nju3Qf8FOioqioifwAuVNU0n3b/APqo6rAT2dO8TbQevU24uvfVdCzqSOHeQvZv2s/+Xfs5UnyEw40Pc7TtUbKaZpEpmeQV5VXYT7RElxKYsqLjOxryHo1iG1X78zMMwwgH6vTISkTigMHAo2WK5gPDK2k2DFjiFSoP84DfACnAZk+d+WXazQOuEZFYVa1g26ij04FYOq4Yyyu5r5DTwD0ipmkMCQkJNI9vTttWbUlrmuYEJr51hVNwrRu1plnDZuZsYBiGUQXCXqyA1kA0sLfM/b3AmZW0aQ/sqKC+t2yz57yggjoxnmfu9i0QkRuBGz2Xed+89f4a3jpeXkghBz3/bWbzSd5SnaI1kBFqI8IE+yyOY5/FceyzOE4PfxtGglh5KTtfKRXcO1n9sverUsfdUH0GeAZARJb7O5Sta9hncRz7LI5jn8Vx7LM4jogs97dtJMxBZQBFuJGQL20pP9rysqeS+vi0qaxOIXDAL0sNwzCMGiHsxUpV84EVwPgyReOBTytptgw4XUQalqm/C9jiU6fsNOJ4YPmJ1qsMwzCM2ifsxcrD48C1InK9iPQSkSeBDsBTACLysIh86FP/JSAbmCMip4jIhcBdwOM++6yeAjqKyGxPn9cD11LekaMingnO26oT2GdxHPssjmOfxXHssziO359FRLiuQ8mm4F/hNgWvAW5X1cWesjnAaFVN8anfF/gLblPwQZw4zapgU/ATHN8U/IeqbAo2DMMwapeIESvDMAyj/hIp04CGYRhGPcbEyjAMwwh7TKyqgYhME5HNIpIrIitE5PRQ2xQKRGSmiGiZY0+o7aoNROQMEXlLRHZ63ve1ZcrF8/nsEpEcEflIRPqEyNwapQqfxZwKviefhcjcGkNE7haRL0UkS0T2i8jbInJKmTr14ntRxc/Cr++FiVUV8Qmm+ztgIM5t/n0R6RRSw0LHOpyzi/foG1pzao0mOAefnwM5FZT/CvgFcBsuiPI+4AMRSag1C2uPk30W4KLE+H5Pzqkd02qV0cBfceHfxuL2ai4QkZY+derL92I0J/8swJ/vharaUYUD+Bz4e5l7G4CHQ21bCD6LmcCaUNsR6gM4Clzrcy24MF33+tyLB44AN4Xa3tr8LDz35gDvhNq2EHwWTXCBDM6z70XpzyKQ74WNrKqATzDdsoFvTxRMt67TxTP9s1lEXhGRLqE2KAxIxUVFKfmeqAumvJj6+z0ZKSL7RGS9iPxdRCpPAld3SMDNWh30XNfn70XZz8JLtb8XJlZV40TBdMuGbKoPfI7bQD0RuAH3GXwqIq1CaVQY4P0u2PfEMRe4GhiHmwIbAiwUkQYhtarmeRJYjYuSA/X7e1H2swA/vxeRFMg2HKhuMN06iaq+73vtWRzdBFyDizZS37HvCaCqr/hcfiMiK4CtwLnAf0JjVc0iIo/jErmOVNWiMsX16ntR2Wfh7/fCRlZVw59guvUGVT0KfAuknaxuHcfrEWnfkwpQ1V241D118nsiIk8AlwFjVXWTT1G9+16c4LMoR1W/FyZWVUD9C6Zbb/AEDO5JmRxg9ZDNuB+mku+J57M5HfueICKtgSTq4PfEE6/0ctyP8/dliuvV9+Ikn0VF9av0vbBpwKrzOPCCiHwBfALcjE8w3fqEiDwKvA1sw/11eD/QGPhXKO2qDUSkCdDNcxkFdBKRAUCmqm4TkdnAvSLyPbAeuA/nKfdSCMytUU70WXiOmcAbuB+hFOBhnMv2m7Vsao0iIn8BrgImAwdFxDuCOqqqR1VV68v34mSfhec7MxN/vhehdm2MpAOYhksxkocbaZ0RaptC9Dm8ggv8mw/s9Hzxeofarlp676Nx6wxljzmecvH8Y9wN5AIfA6eE2u7a/ixwrtnzPD9C+bg1iTlAcqjtroHPoaLPQIGZPnXqxffiZJ9FIN8LC2RrGIZhhD22ZmUYhmGEPSZWhmEYRthjYmUYhmGEPSZWhmEYRthjYmUYhmGEPSZWhmEYRthjYmUYQaKChHIVHVtEJKWiZIUhsjmljH2jfco+EpGlVejjPp/2O2rSXqP+YhEsDCN4DCtz/SbwFW4zqJc83MbQYcAPtWNWlfgt8C7wnR9tn8Ml07sf6B9MowzDi4mVYQQJVS2VmltE8oCMsvc9hFt69x8qsfOkqOpOYKeI7A+yTYZRgk0DGkYtU9E0oIjMEZEdIpIuIp+KSI6IrBORcz3ld3imELNE5H8i0qZMnzEicreIfC8ieSKyS0Qe8wRMDdTeM0VkpYhki8gaEZkcaJ+GUV1MrAwjfGgKPA/8A7gAFz/tDRF5DBgD3AJM97z+S5m2/8YFR30JlxfoYWAq8GKANnXFJdB7HLgQN4X5uoh0O2ErwwgyNg1oGOFDAnCzqi4GEJFduDWvSbhAwUWe+6cAt4lItKoWicjpwE+Aa1T1eU9fC0QkE/i3iAxQ1dV+2tQaF7B5g+fZK3GCdQnwOz/7NIxqYyMrwwgfjnmFyoM3F9ACLZ119nvcH5qJnusJuAjWb3imA2NEJAaY7yk/IwCbNniFCkBV9+FGfJ0C6NMwqo2NrAwjfDjke6Gq+SICcLBMvXzP2bse1RaIw+VHqohWAdiUWcG9PJ9nG0atYGJlGJHPAVyOpNMrKd9Vi7YYRo1gYmUYkc9cYAbQTFU/DLUxhlETmFgZRoSjqh+JyMs4L73HgS+AYlzK8HOAGaq6PoQmGkbAmFgZRt3gSuA24DrgXty60hZcCvG9oTPLMIKDpbU3jHqMiKQAm3F7sp4HirSaPwrivECigWeBcaraMdh2Goa5rhuGAU5oCoBRfrS919P26qBaZBg+2MjKMOoxIhIH9PO5tU5Vj1Szj0QgyXOZr6pfB8s+w/BiYmUYhmGEPTYNaBiGYYQ9JlaGYRhG2GNiZRiGYYQ9JlaGYRhG2GNiZRiGYYQ9/x8WQsBb/zllXAAAAABJRU5ErkJggg==\n",
      "text/plain": [
       "<Figure size 432x288 with 1 Axes>"
      ]
     },
     "metadata": {
      "needs_background": "light"
     },
     "output_type": "display_data"
    }
   ],
   "source": [
    "plt.figure(figsize=[6,4])\n",
    "plt.plot(x,solTchem[:,Header.index('OH')],'m',label='OH')\n",
    "plt.plot(x,solTchem[:,Header.index('O')],'r',label='O')\n",
    "plt.plot(x,solTchem[:,Header.index('O3')],'k',label='$\\mathrm{O_3}$')\n",
    "# plt.plot(x,solTchem[:,Header.index('NO')])\n",
    "plt.plot(x,solTchem[:,Header.index('NO2')],'g',label='$\\mathrm{NO_2}$')\n",
    "plt.xlim([0,25])\n",
    "plt.ylim([0,2.e-6])\n",
    "plt.ylabel('Concentration '+ r'[$\\mathrm{\\frac{mol}{m^{-3}}}$]',fontsize=16)\n",
    "plt.xlabel('Time [h]',fontsize=16)\n",
    "plt.xticks(fontsize=14)\n",
    "plt.yticks(fontsize=14)\n",
    "plt.legend()\n",
    "plt.savefig('species_concentration.pdf',bbox_inches='tight')"
   ]
  },
  {
   "cell_type": "code",
   "execution_count": null,
   "metadata": {},
   "outputs": [],
   "source": []
  }
 ],
 "metadata": {
  "kernelspec": {
   "display_name": "Python 3 (ipykernel)",
   "language": "python",
   "name": "python3"
  },
  "language_info": {
   "codemirror_mode": {
    "name": "ipython",
    "version": 3
   },
   "file_extension": ".py",
   "mimetype": "text/x-python",
   "name": "python",
   "nbconvert_exporter": "python",
   "pygments_lexer": "ipython3",
   "version": "3.8.11"
  }
 },
 "nbformat": 4,
 "nbformat_minor": 4
}
